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Abstract 

We report the trigonometric parallax of IRAS 07427-2400 with VERA to be 0.185 ± 

0.027 mas, correspond.ing to a distance of 5.4llg gg kpc. The result is consistent with the 
previous result of 5.32fQ'42 kpc obtained by Choi et al. (2014) within error. To remove 
the effect of internal maser motions (e.g., random motions), we observed six maser features 
associated with IRAS 07427-2400 and determined systematic proper motions of the source 
by averaging proper motions of the six maser features. The obtained proper motions are 
ijtaCosS, pLs) = (—1.79 ± 0.32, 2.60 ± 0.17) mas yr“^ in equatorial coordinates, while Choi 
et al. (2014) showed (fj,aCosS, fis) = (—2.43 ± 0.02, 2.49 ± 0.09) mas yr“^ with one maser 
feature. Our astrometry results place the source in the Perseus arm, the nearest main arm in 
the Milky Way. Using our result with previous astrometry results obtained from observations 
of the Perseus arm, we conducted direct (quantitative) comparisons between 27 astrometry 
results and an analytic gas dynamics model based on the density-wave theory and obtained 
two results. First is the pitch angle of the Perseus arm determined by VLBI astrometry, 11.1 
±1.4 deg, differing from what is determined by the spiral potential model (probably traced 
by stars), ^ 20 deg. The second is an offset between a dense gas region and the bottom 
of the spiral potential model. The dense gas region traced by VLBI astrometry is located 
downstream of the spiral potential model, which was previously confirmed in the nearby 
grand-design spiral galaxy M51 in Egusa et al. (2011). 
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1. Introduction 

In 1964, Lin & Shu (1964) proposed the density- 
wave theory as a way of overcoming the winding 
dilemma, in which a spiral arm in a disk galaxy is 
destroyed within a time scale of several galactic ro¬ 
tations due to the differential rotation. Since then 
the density-wave theory has been used to discuss 
the evolution and dynamics of a spiral arm in a disk 
galaxy. For instance, Fujimoto (1968) and Roberts 

* Last update: Januaiy 19, 2007 


(1969) examined non-linear gas motions perturbed 
by the density wave (spiral arm) numerically, and 
found the “Galactic shock” where a velocity jump 
(non-circular motion) occurred and the gas accu¬ 
mulated upstream of the spiral potential (see fig. 5 
in Roberts 1969). Burton (1973) and Mel’nik et al. 
(1999) schematically showed the systematic non¬ 
circular motion based on the density-wave theory. 

In contrast to the density-wave theory, another 
dynamics model, the recurrent transient spiral, has 
been proposed based on N-body simulations with 
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or without a gas disk (e.g.. Miller et al. 1970; 
Hohl 1971; Hockney & Brownrigg 1974; James & 
Sellwood 1978; Sellwood & Carlberg 1984; Baba 
et al. 2010; Wada et al. 2011; Fuji! et al. 2011; 
Grand et al. 2012; Baba et al. 2013). In the recur¬ 
rent transient spiral, the spiral arm is destroyed and 
regenerated over several rotational periods. Wada 
et al. (2011) reported that gas flows converge near 
the bottom of the spiral potential with random mo¬ 
tions, which is different from the result based on 
the density-wave theory. 

Ever since the 1960s, the two models have been 
proposed and suggested by simulations and obser¬ 
vations. However, the previous observations were 
mainly based on line-of-sight velocities (ID infor¬ 
mation) and apparent 2D positions on the celestial 
sphere, which were not enough to conduct quan¬ 
titative comparisons between the observations and 
the simulations in terms of both velocity and spatial 
resolutions. 

Since the 2000s the Very Long Baseline 
Interferometry (VLBI) technique has allowed us 
to routinely conduct Galactic-scale astrometry and 
to obtain accurate information about 3D positions 
and velocities of Galactic masers in the Milky Way 
(e.g., Reid & Honma 2014a; Reid et al. 2014b). At 
the moment, we can conduct direct (quantitative) 
comparisons between the astrometry results and the 
dynamical models for understanding the dynamics 
and evolution of the spiral arm, as well as those of 
the Milky Way. 

VLBI astrometry observations have revealed 
systematic peculiar (non-circular) motions in the 
Perseus arm (e.g., Sakai et al. 2012, 2013; Choi 
et al. 2014). The Perseus arm traced by stars and 
gas is one of the Milky Way’s main arms (e.g., 
Churchwell et al. 2009). Sakai et al. (2012) 
and (2013) studied the 3D structure and kinematics 
of the Perseus arm using VLBI astrometry results. 
The sources in the Perseus arm were moving sys¬ 
tematically toward the Galactic Center and lagged 
behind the Galactic rotation. 

Sakai et al. (2012) determined averaged pecu¬ 
liar motions of (f/mean, ^mean) = (11 ± 3, -17 ± 
3) km s“^ with seven sources in the Perseus arm 
with ^ 4 -(t significance. Note that U and V are di¬ 
rected to the Galactic center and the Galactic rota¬ 
tion, respectively. The systematic peculiar motions 
were also confirmed by Choi et al. (2014), who 
determined averaged peculiar motions of (C/mean, 
Knean) = (9.2 ± 1.2, -8.0 ± 1.3) km s-i with 25 


sources in the Perseus arm. The systematic pecu¬ 
liar motions are consistent with the Galactic shock 
proposed by Lujimoto (1968) and Roberts (1969). 
Thus, the Perseus arm is a good target for testing 
the density-wave theory. 

IRAS 07427-2400, a high-mass star-forming re¬ 
gion (Qiu et al. 2009), had the largest Galactic lon¬ 
gitude {I ~ 240°) in the Perseus arm in Choi et al. 
(2014). The source played an important role in the 
study of the large-scale structure and kinematics of 
the Perseus arm (e.g., determination of the pitch an¬ 
gle). However, Choi et al. (2014) observed only 
one maser feature to discuss the systemic motion 
of the source (e.g.. Galactic rotation), while Qiu 
et al. (2009) confirmed the bipolar molecular out¬ 
flow of ^^CO(J=2-l) around the source. To discuss 
the systemic motion of the source precisely, more 
maser features should be observed to remove the 
effect of internal maser motions (e.g., outflow mo¬ 
tions). 

In addition, a different distance has rarely been 
measured for the same source with independent 
VLBI arrays (e.g., d = 5.03 ± 0.19 kpc for 
G048.60 h- 0.02 in Nagayama et al. 2011; d = 
10-75lo;55 kpc for G048.60 h- 0.02 in Zhang et 
al. 2013), although the same distances have 
been measured for the same sources in most 
cases. Therefore, independent distance measure¬ 
ments with independent arrays are important to 
confirm astrometry results precisely. 

In this paper, we will report astrometry results of 
IRAS 07427-2400 with VERA for discussing the 
large-scale structure and kinematics of the Perseus 
arm more precisely. Using our results with previ¬ 
ous astrometry results, we will compare the results 
with an analytic gas dynamics model for testing the 
density-wave theory. 

2. Observations and Data Reduction 

2.7. VLBI Observations with VERA 

We carried out VERA astrometry observations 
of H 2 O maser at a rest frequency of 22.235080 
GHz to measure the trigonometric parallax and 
proper motions of IRAS 07427-2400. Note that 
IRAS 07427-2400 was observed as a part of the 
Outer Rotation Curve project with VERA (as ex¬ 
plained in Sakai et al. 2012), aiming to understand 
the mass distribution of the Milky Way especially 
for the outer disk based on an accurate rotation 
curve. The observations procedure has been de- 
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No. ] 


Table 1. Observations. 


Epoch 

Date 

Time Range 
(UTC) 

Beam 

(mas) 

Image r.m.s. Detected maser feature 
(mjy/beam) 

A 

Jan 18, 2012 

11:00-18:05 

1.7x0.8@157° 

124 

1,2, 3, 4 

B 

Feb 28, 2012 

08:20-15:25 

1.7x0.8@157° 

180 

1,2, 3, 4 

C 

Apr 24, 2012 

04:40-11:45 

1.7x0.8@162° 

154 

1, 2, 3, 4, 5 

D 

Jun 17, 2012 

01:10-08:15 

1.7x0.7@157° 

356 


E 

Sep 10, 2012 

19:30-02:35 

1.9x0.9@162° 

384 


F 

Jan 27, 2013 

09:55-18:20 

1.9x0.8@162° 

69 

1, 2, 3, 4, 5, 6 

G 

Mar 20, 2013 

06:30-14:55 

1.8x0.9@158° 

159 

1, 2, 3, 4, 5, 6 

H 

Jun 02, 2013 

01:35-10:00 

1.9x0.8@159° 

249 

1,2, 3, 4, 5,6 

I 

Sep 16, 2013 

18:35-03:00 

1.9x0.8@160° 

103 

1,2, 3, 4, 5,6 


Table 2. Source Data. 


Source Name 

R.A. 

(J2000.0) 

Deck 

(J2000.0) 

S.A.* 

(°) 

Flux Density 

(Jy) 

Note 

IRAS 07427-2400 

07^44“5I?9205 

_24^07'4I"457 


II - 

-27 

H 2 O masers 

J0745-245I 

07‘'45“I0f263445t 

-24'^5I'43"7668I'f 

0.74 

0.07- 

-0.24 

Phase-reference source 


'Separation angle between the maser and reference sources. 
^The positions are based on Petrov et al. (2011). 


scribed in previous VERA astrometry papers (e.g., 
Sakai et al. 2012). 

In table 1, the observation periods and synthe¬ 
sized beam sizes are listed with image noise lev¬ 
els and information about detected maser features. 
In table 2, maser and reference source information 
are listed for phase-referencing observations (e.g., 
tracking-center positions, separation angle between 
the two sources, and flux densities). Note that we 
observed one or two bright continuum source(s) 
in each observation to conduct a clock offset cal¬ 
ibration of each VERA station. The observed con¬ 
tinuum sources were 3C84, DA193, M87, OJ287, 
and BLLAC from which we selected one or two 
source(s) in each observation. 

As for an on-source time of the maser source in 
each observation, ^ 2 hours were assigned in the 
total observation time of ^ 7 hours in the five for¬ 
mer observations, while ~ 3 hours were assigned 
in the total observation time of ^ 8 hours in the 
four later observations (as shown in table 1). Note 
that we observed two pairs of maser and reference 
sources for VLBI astrometry observations. One of 
them is listed in table 2 and reported in this paper. 
The other pair, IRAS 07207-2400 and J0729-1320, 
will be reported in another paper. Eor velocity reso¬ 
lution of the maser source we set 0.42 km s“^ with 
a velocity coverage of ~ 215 km s“^ , instead of a 


high resolution mode of 0.21 km s ^ with a veloc¬ 
ity coverage of ^ 108 km s“^ (e.g., explained in 
Sakai et al. 2014). 

2.2. Data Reduction 

We conducted general phase-referencing anal¬ 
ysis with AIPS (Astronomical Image Processing 
System, NRAO) in the same manner as Kurayama 
et al. (2011). The main analysis procedures are de¬ 
scribed as follows (see more detail in the Appendix 
of Kurayama et al. 2011): 

1. Amplitude calibration for line (maser) 
and continuum (phase-reference and clock- 
calibration) sources 

2. Accurate recalculation of a tracking model 
used in the Mitaka EX correlator 

3. Binding of frequency channels for the con¬ 
tinuum sources 

4. Calibration of clock parameters with the 
clock calibrator 

5. Calibration of clock parameters with the ref¬ 
erence source 

6. Imaging of the reference source after the self 
calibration 

7. Phase referencing for the maser source (sub¬ 
tracting phases of the reference source from 
those of the maser source) 
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8. Dual-beam calibration with the “horn-on- 
dish” method (Honma et al. 2008) 

9. Calibration of the Doppler effect for the 
maser source 

10. Making a CLEANed image of the maser 
source with the phase referencing in each ve¬ 
locity channel 

11. Determining an absolute position and a flux 
density of the maser source in each velocity 
channel. 

Step 2 (accurate recalculation of a tracking model) 
was needed, since the tracking model used in the 
correlator did not have sufficient accuracy for as¬ 
trometry. For step 11 (determining an absolute po¬ 
sition and a flux density of the maser), we con¬ 
ducted an elliptical Gaussians fit with the task “jm- 
flt” in AIPS. 

We regarded a maser as detected if it achieved a 
high SNR (Signal-to-Noise Ratio) of five or more. 
We traced a maser detected in the same velocity 
channel and its position movements during the ob¬ 
servations to determine the parallax and proper mo¬ 
tions. We determined the proper motions using 
masers detected with three epochs or more. On 
the other hand, we determined the parallaxes using 
masers detected in at least two continuous velocity- 
channels with five epochs or more. 

On determinations of the parallax and proper 
motions, we used the “VERA Parallax”, one of the 
tasks of “VEDA (VEra Data Analyzer)”, a data an¬ 
alyzing software developed at NAOJ. In the task, 
we assumed that source motions can be modeled 
by a combination of linear proper and sinusoidal 
parallax motions. 

3. Results 

3.1. Maser spatial/velocity distributions in IRAS 
07427-2400 

During the observations over a period of one 
year and eight months, we detected maser emis¬ 
sions except for epochs D and E as shown in table 
1. For epochs D and E we could not make good 
quality images due to high image noise levels (e.g., 
table 1). Figure 1 shows a typical maser spectrum 
taken at epoch F. 

Figure 2a shows a large-scale view of IRAS 
07427-2400 with which the bipolar molecular out¬ 
flow of ^^CO(J=2-l) is associated (Qiu et al. 2009). 
Figure 2b represents the magnified view of fig. 2a, 
and there are three continuum (1.3 mm) peaks as 
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Fig. 1. Scalar-averaged cross-power spectium of IRAS 

07427-2400 taken at epoch F with several baselines av¬ 
eraged. 

MMl, MM2, and MM3 (Qiu et al. 2009) in fig. 
2b. Figures 2c and 2d show VERA observation re¬ 
sults in the same area of the rectangular area in fig. 
2b. Fabeled numbers in figures 2c and 2d repre¬ 
sent maser features listed in table 1. The maser 
feature is recognized as a maser cluster, which 
consists of closest maser spots in continuous ve¬ 
locity channels. In figure 2c, one can see maser 
spatial/velocity distributions and observed motions 
relative to the reference source listed in table 2. 

3.2. Trigonometric parallax of IRAS 07427-2400 

As a result of the data analyses, we detected 
three maser features to determine the parallaxes as 
shown in table 3. Using the three maser features 
consisting of 10 spots, we determined the parallax 
of IRAS 07427-2400 to be 0.185± 0.015 mas in 
table 3. The position errors shown in table 3 for 
right ascension (R. A.) and declination (Deck) were 
introduced so that the reduced yf became unity, 
since systematic error is the dominant source of er¬ 
ror in VFBI astrometry observation (e.g., Sanna et 
al. 2012). The position error of R.A., 0.069 mas, 
is smaller than that of Deck, 0.204 mas, which 
is consistent with previous VERA results. As for 
the obtained parallax error (0.015 mas), it may be 
underestimated since the 10 spots may not be in¬ 
dependent of each other. For a conservative par¬ 
allax estimation, we multiplied the obtained error 
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Fig. 2. (a) Large-scale view of IRAS 07427-2400 with which the bipolar outflow of ^®CO(J=2-l) is associated (Qiu et 
al. 2009). Red and blue contours represent redshifted and blueshifted ^®C0(J=2-1) lobes, respectively. Redshifted and 
blueshifted lobes were integrated from 74 to 78 km s“^ and from 58 to 62 km s“^, respectively, (b) Magnified view of 
the rectangle of fig. 2a. Contours show continuum (1.3 mm) emissions in which there are three 1.3-mm peaks as MMl, 
MM2, and MM3 (Qiu et al. 2009). (c) Maser spatial and velocity distributions (as colored circles) with observed motions 
(as green aiTows) relative to the phase reference source in table 2. The distributions were made based on epoch F in table 
1. We adopted our distance measurement of 5.41 kpc to show the observed motions in absolute velocities. The numbers 
labeled correspond to the maser features listed in table 1. Nominal origin is set to the maser detection position in table 2. 
(d) Same as (c), but with systematic proper motions subtracted from fig. 2(c) (i.e., internal motions of masers) (see text). 
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Table 3. Pai-allax fits.* 


Feature 

Flsr 
km s“^ 

Nepoc/is 

Epochs 

Parallax(Error) 

(mas) 

Errors 

R.A. Dec. 

(mas) 

1 

68.7 

6 

ABC-h-hFG-hI 

0.171(0.039) 

0.049 

0.162 


68.2 

7 

ABC-l-l-FGHI 

0.162(0.047) 

0.066 

0.200 


67.8 

7 

ABC-l-l-FGHI 

0.184(0.064) 

0.094 

0.171 


67.4 

7 

ABC-l-l-FGHI 

0.182(0.073) 

0.110 

0.180 


67.0 

6 

A-hC-h-hFGHI 

0.181(0.075) 

0.109 

0.195 

2 

66.6 

7 

ABC-l-l-FGHI 

0.193(0.029) 

0.041 

0.230 


66.1 

7 

ABC-l-l-FGHI 

0.205(0.020) 

0.028 

0.201 


65.7 

6 

A-hC-h-hFGHI 

0.163(0.035) 

0.048 

0.217 

3 

64.9 

7 

ABC-l-l-FGHI 

0.225(0.018) 

0.025 

0.296 


64.4 

6 

ABC-h-hFG-hI 

0.218(0.030) 

0.037 

0.132 

Combined fit for 10 spots 



0.185(0.015) 

0.069 

0.204 

Final 




0.185(0.027) 




* A combined fit was done to a data set of 10 spots (see text). The error of the final parallax was estimated to multiply 
the error of the combined fit by where A^spot and Wfeature are the numbers of maser spots and features, 

respectively (see also text). 


Table 4. Determination of systematic proper motions for IRAS 07427-2400. 


Feature 

Vlsr 
km s“^ 

Proper Motions* (Error) 
PoCOSl) fis 

(mas yr“^) (mas yr“^) 

Note 

1 

69.9 ~ 66.6 

-2.55 

2.51 


2 

66.6 ~ 65.7 

-1.45 

2.10 


3 

64.9 ~ 64.4 

-0.53 

3.31 


4 

67.0 ~ 66.6 

-2.44 

2.83 


5 

70.4 ~ 68.7 

-1.50 

2.32 


6 

63.6 ~ 62.3 

-2.26 

2.51 


Mean 

66.4^ 

-1.79 (0.32)» 

2.60 (0.17)f 

Systematic proper motions 


* Proper motions were determined by adapting a parallax of 0.185 mas. 
^ The value was determined by averaging the results of the six features. 
^ The standai'd en'ors. 


by '^here A^spot and iVfeature are the 

numbers of maser spots (-/Vspot = 10) and features 
(A^feature = 3), respectively. 

We obtained the hnal parallax to be 0.185 
± 0.027 mas, corresponding to a distance of 
5.41 ^Qgg kpc. The distance is consistent with the 
previous parallactic distance of 5.321 q 42 kpc de¬ 
termined by Choi et al. (2014) within error. On the 
other hand, the distance is smaller than the previ¬ 
ously estimated source distance of 6.9 kpc, based 
on the kinematic distance derived from ^^CO(J=l- 
0) observations (Wouterloot & Brand 1989). The 
difference between the distances can be explained 


using peculiar (non-circular) motions of the source. 
We will further discuss the peculiar motions in 
chapter 4. 

Figure 3 shows combined ht results in which 
maser motions are modeled by a combination of 
common sinusoidal parallax and discrete linear 
proper motions. Figures 3a and 3b show the com¬ 
bined ht results of each maser spot with the h- 
nal parallax result, hxed in directions of R.A. and 
Deck, respectively. Error bars in each panel of hg- 
ure 3 represent the systematic position errors ex¬ 
plained previously. Figures 3c and 3d show the 
sinusoidal parallax motions after subtracting the 
proper motions from hgures 3a and 3b in direc- 
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Fig. 3. Maser positional evolutions and combined-fit results for the 10 maser spots listed in table 3 (see text). Error bars 
represent position errors resulting from systematic errors, which are introduced so that the reduced becomes unity. 
Different colors coiTespond to respective features, (a) The maser’s positional evolutions in right ascension as shown by 
circles. Solid lines show the parallax and proper motion fits (see text), while dashed lines represent only the proper motion 
fits, (b) Same as (a), but in declination, (c) Parallax motions (as the circles) in right ascension after subtracting the proper 
motions from (a). The solid curve shows the modeled parallax motion, (d) Same as (c), but in declination. 
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tions of R.A. and Deck, respectively. Large devia¬ 
tions from the model are seen in Deck (figure 3d), 
which have been shown in previous VLBI astrom¬ 
etry results (e.g., Sakai et al. 2012). The main rea¬ 
son for the deviation is thought to be tropospheric 
zenith delay residuals as discussed in Honma et al. 
(2007). 

3.3. Systematic proper motions of IRAS 07427- 
2400 

Since observed proper motions include systemic 
(e.g., Galactic rotation) and internal maser (e.g., 
outflow) motions, we should separate both motions 
for discussing the systemic motion of IRAS 07427- 
2400. Table 4 displays observed proper motions of 
several maser features that were detected with three 
epochs or more. Note that the proper motion com¬ 
ponents were determined by adapting a common 
parallax of 0.185 mas. 

In fig. 2c we can see clear northwest motions, 
which represent the Galactic rotation. The velocity 
range of maser features in fig. 2c is not consistent 
with velocity ranges of redshifted and blueshifted 
lobes of molecular outflow (Qiu et al. 2009). 
Thus, we assumed that internal maser motions can 
be regarded as random motions, which means that 
we can determine the systematic proper motions by 
averaging proper motions of several maser features. 
Averaging of multiple features is essential for ob¬ 
taining the systematic proper motions, since table 4 
shows that a difference of observed proper motions 
is up to ^ 2 mas, corresponding to ^ 50 km s“^ at 
a distance of 5.4 kpc. The difference is caused by 
internal maser motions. 

As a result of the data averaging in table 4, we 
determined the systematic proper motions to be 
{fiaCosS, p,s) = (—1.79 ± 0.32, 2.60 ± 0.17) mas 
yr“^ in equatorial coordinates. Note that the ob¬ 
tained errors were determined by dividing standard 
deviations by a factor of where n is the number 
of measurements (e.g., n = 6 applied in this case as 
shown in table 4). Table 4 shows a mean veloc¬ 
ity (Llsr) of ~ 66 km s“t^ which is ~ 2 km s“t 
smaller than a velocity of CS(J=2-1) emission in 
IRAS 07427-2400 (Bronfmann et al. 1996). The 
difference indicates that the systematic proper mo¬ 
tions include errors of at least ^ 2 km s“^. In fact, 
the obtained errors in table 4, 0.32 mas and 0.17 
mas, are converted to ^ 8 km s“t and ^ 4 km s“t 
at a distance of 5.4 kpc, respectively. 

Figure 2d represents internal maser motions af¬ 


ter subtracting the systematic proper motions from 
fig. 2c, which in fact seems to be the random mo¬ 
tions. We will discuss the systemic motion of IRAS 
07427-2400 using the resultant systematic proper 
motions in chapter 4. 

3.4. Peculiar motions of IRAS 07427-2400 

As the next step, we determine the peculiar (non¬ 
circular) motions of IRAS 07427-2400. Full-space 
(3D) motions of the source can be determined from 
the observed 3D position, systematic proper mo¬ 
tions, and systemic velocity of the source. For the 
systemic velocity we refer to 68.0 ± 5.0 km s“^ ob¬ 
tained from CS(J=2-1) observations (Bronfmann et 
al. 1996). Calculation procedures for determining 
the 3D motion are itemized as follows: 

1. Velocity conversion from the LSR velocity 
(I^Sr) to the heliocentric velocity (wheiio) 

2. Coordinates conversion from equatorial co¬ 
ordinates to Galactic coordinates 

3. Conversion from the proper (apparent) mo¬ 
tions to absolute motions using the parallax 
result (tt) 

4. Corrections of the peculiar solar motions, 
the Galactic constants (i?o and ©o), and the 
Galactic rotation curve (0(i?)) 

Note that detailed explanations for the above 
steps are summarized in the Appendix 1. 

In Step 1 (Velocity conversion) Vlsr = 68.0 ± 
5.0 km s“^ is converted to Uheiio = 86.4 ± 5.0 km 
s“^ for the source. For Steps 2 and 3 we use the 
matrix multiplication shown in the Appendix 1. In 
Step 4 we adopt peculiar solar motions of (Uq, Vq, 
Wq) = (11.1 ± 1, 12.24 ± 2, 7.25 ± 0.5) km s"! 
(Schonrich, Binney, & Dehnen 2010) and Galactic 
constants of {Rq, ©q) = (8.33 kpc, 240 km s“^) 
referred from Gillessen et al. (2009) and Reid & 
Brunthaler (2004). For the Galactic rotation curve 
we simply assume a flat rotation curve as ©(i?) = 
© 0 . 

As a result, we obtain the peculiar motions of 
IRAS 07427-2400 to be {U, V, W) = (-8 ± 8, 
0±9, 1 ± 11) km s“^. Directions of the pecu¬ 
liar motions are toward the Galactic center (U), the 
Galactic rotation (V), and the north Galactic pole 
(W). Note that here the obtained errors are esti¬ 
mated considering contributions from the parallax 
(St!-), the proper motions (Sfj,^cosS, sys¬ 

temic velocity (<),;heiio) basically the same man¬ 
ner as Johnson & Soderblom (1987). A detailed 
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Fig. 4. (a) Schematic face-on view of the Milky Way (Taylor Sz Cordes 1993 with updates) on which our astrometry result of IRAS 
07427-2400 and previous astrometry results are superimposed. Rq = 8.33 kpc (Gillessen et al. 2009) was assumed, and the 
Sun is located at the origin. Red and black circles represent the sources in the Perseus and Local arms, respectively, while 
green ones, refeiTed from Sakai et al. (2012), show the sources in the other regions, (b) Galactocentric azimuth (/3) versus 
Galactocentric distance (R) diagram on which the astrometiy results are superimposed. Note that R axis is shown in the log 
scale. The solid line based on a logarithmic spiral model (see eq. 1 in text) was fitted to the sources in the Perseus arm, while 
the dashed line shows the previous result obtained by Sakai et al. (2012). 


explanation of the error calculation is described in 

the Appendix 1. 

4. Discussion 

4.1. Large-scale structure of the Perseus arm 

In this section, we determine a pitch angle and 
a reference position of the Perseus arm based on 
our parallax result and previous VLBI astrometry 
results. First, a combination of coordinates (e.g.. 
Galactic longitude and latitude) and the parallax of 
IRAS 07427-2400 tells us the absolute source po¬ 
sition on the Galactic plane as shown in fig. 4a. 
Note that we assumed a distance to the Galactic 
center (i.e., the Galactic constant Rf) of 8.33 kpc 
referred from Gillessen et al. (2009). The source 
is positioned only ~ 6 pc away from the disk to¬ 
ward the north Galactic pole. Second, based on the 
source location and the longitude-velocity {l-v) di¬ 
agram of CO(J=1-0) emission (Dame et al. 2001), 
the source is positioned in the Perseus arm (e.g., 
figures 4a and 4b, and fig. 3 in Dame et al. 2001). 


A shape of the spiral arm is often described with a 
logarithmic-spiral model (e.g., Reid et al. 2009a) 
as 

ln(ii!/i?ret) =-(/3-/3ref)tam, (1) 

where i is the pitch angle of the spiral arm and /3 
is Galactocentric azimuth (/3 is defined as 0 deg 
toward the Sun from the Galactic Center and in¬ 
creasing clockwise). (i?ret, /3ref) are reference po¬ 
sitions of the spiral arm. Since {Rref, /3ref) can 
have an infinity of combinations in the log-spiral 
fit, we regard /3ref as 0 deg in this paper for sim¬ 
plicity. Finally, we apply the log-spiral model to 
our result and the previous ones (as summarized in 
fig. 4b and table 5), then we can determine a pitch 
angle and a reference position of the Perseus arm. 

As a result, we obtain the parameters, {i, Rref), 
to be (11.2 ± 1.4 deg, 10.4 ±0.1 kpc) using 
an unweighted least-squares fit toward 27 sources 
in the Perseus arm. If we exclude W49N and 
G48.60 h- 0.02 located in the Perseus arm of the first 
Galactic quadrant from the fit, (i, i?ret) = (12.2 ± 
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2.1 deg, 10.4 ±0.1 kpc) are obtained, which are 
consistent with the results obtained for 27 sources 
within errors. The former pitch angle is also con¬ 
sistent with the previous result of 9.9 ±1.5 deg 
determined by Choi et al. (2014) using 25 sources 
in the Perseus arm within error. 

On the other hand, the pitch angle (11.2 ± 1.4 
deg) is not consistent with the previous result of 
17.8 ±1.7 deg determined by Sakai et al. (2012) 
using seven sources in the Perseus arm. The pitch 
angle difference (11.2 ± 1.4 deg or 17.8 ± 1.7 deg) 
may originate in the difference of the source num¬ 
bers (27 or 7). However, here we cannot choose 
one as the best solution, since an observation gap 
is seen in the Perseus arm (e.g., 28.8 deg < /3 < 
81.4 deg in fig. 4b). More astrometry observations 
especially for the gap region would help us to judge 
whether a spur or bifurcation of the spiral arm, as 
seen in external disk galaxies, exists or not in the 
Perseus arm. 
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No. ] 


Table 5. Peculiar motions in the Perseus and Local arms.* 


Source 

1 

b 

R 

/3+ 

u 

V 

W 

Arm 

Ref.i* 


deg 

deg 

kpc 

deg 

km s“^ 

km 

km 



EC 95 

31.56 

5.33 

7.98 

1.6 

-1±4 

2i3 

2i0.5 

Local 

6 

W49N 

43.17 

0.01 

/■OW_o,44 


-17±14 

-24il3 

-5il4 

Perseus 

32 

G48.60h-0.02 

48.61 

0.02 

o 1 ^+0.41 
^■^^-0.34 

8i.4+:l 

-5±12 

12ill 

6il2 

Perseus 

32 

G59.78h-0.06 

59.78 

0.06 

7 4Q+0-02 

14.4t!];? 

1±2 

-li3 

-4il 

Local 

29 

ONI 

69.54 

-0.98 

7 QO+0.01 

/■oZ_o 01 


7±5 

-6i5 

5i5 

Local 

17,21,30 

G74.03-1.71 

74.03 

-1.71 

8.04 

11.0t°:| 

8±2 

-6i5 

10i2 

Local 

30 

G75.76h-0.33 

75.76 

0.34 

8-20;H4 

24 . 5 + 1 ° 

0±6 

-16i9 

6i5 

Local 

30 

ON2N 

75.78 

0.34 


25 . 7 tli 

-3±6 

-5i5 

li5 

Local 

2,30 

G76.38-0.61 

76.38 

-0.62 

o 1 o+O.Ol 
^■^^-0.01 

8-9+:^ 

3±18 

-lli5 

12il9 

Local 

30 

IRAS 20126H-4104 

78.12 

3.63 

8.15 

1 1 O+0.4 
^ ^■■^-0.4 

6±1 

-13i5 

15il 

Local 

16 

AFGL 2591 

78.89 

0.71 

O 

o--^J_0.02 

23.0+°-? 

-13±4 

-lli7 

-22i5 

Local 

22 

IRAS 20290H-4052 

79.74 

0.99 

8.20 

9.4+;? 

2±3 

-10i5 

6i3 

Local 

22 

G79.87h-1.17 

79.88 

1.18 

8.20 

11.2t+5 

9±10 

-llilO 

3il0 

Local 

30 

NML Cyg 

80.80 

-1.92 

Q '^o+U.UUb 
^■^■^-0.003 

^ ^-+-0.8 

-2±4 

-9i3 

-5i4 

Local 

31 

DR 20 

80.86 

0.38 

8.23 

10 

7±2 

-9i5 

5i2 

Local 

22 

DR 21 

81.75 

0.59 

8.25 

10 4 + 0 ® 

0±2 

-8i3 

7i2 

Local 

22 

W75N 

81.87 

0.78 

8.25 

0 qH-0.5 

Oil 

2i3 

lil 

Local 

22 

G90.2lH-2.32 

90.21 

2.32 

8.36 

+•^-0.1 

-4il0 

-6i5 

6 il 0 

Local 

30 

G92.67h-3.07 

92.67 

3.07 

8.56;Hi 

11 . 0 +:^ 

-16i3 

- 6 il 0 

- 2 i 2 

Local 

30 

AFGL 2789 

94.60 

-1.80 

9 

^■^^-0.11 

19.6t+ 

5i5 

-28i3 

-lli 6 

Perseus 

19 

G094.60-1.79 

94.60 

-1.80 

^■■^^-0.16 

24.5t+ 

12i7 

- 20 i 6 

-15i8 

Perseus 

5 

G095.29-0.93 

95.30 

-0.94 

lU.UZ_o 09 

78 8+0-8 

^O-O_ 0.8 

-9i5 

-7i5 

2i5 

Perseus 

5 

GlOO.37-3.57 

100.38 

-3.58 

' -0.09 

20 . 8 t+ 

13il0 

-3il0 

3il0 

Perseus 

5 

G105.4lH-9.87 

105.42 

9.88 

Q /cr\“l“ 0.02 

8.60I002 

c /;+0.3 

•^•^-0.3 

8 il 0 

-li 6 

-13il0 

Local 

30 

IRAS 22198H-6336 

107.30 

5.64 

0 ^Q + 0.01 

+■^-0.2 

5i2 

-7i5 

llil 

Local 

9,11 

L 1206 

108.18 

5.52 

8-60Io.o3 

4 q +0.6 
+•^-0.5 

0i3 

-8i3 

0 i 6 

Local 

21 

G108.20h-0.58 

108.21 

0.59 

10 S7+0-59 

1WO/_0.40 

23.3t+ 

-llilO 

-15i9 

9ill 

Perseus 

5 

IRAS 22480H-6002 

108.43 

0.89 

^■^^-0.12 

14.6+:« 

12i3 

-23i2 

0i3 

Perseus 

12 

G108.47-2.81 

108.47 

-2.82 

9 84+9.06 

V.O^_0,06 

18.1+:? 

20i7 

-13i6 

-9i7 

Perseus 

5 

G108.59h-0.49 

108.59 

0.49 

Q 41 +0.12 
^■+^-0.10 

14 4+11 
1+-+-1.0 

46i6 

-5i6 

0 i 6 

Perseus 

5 

Cep A 

109.87 

2.11 

8.59;H2 

4 4H-O.3 

2i4 

-2i5 

-5i2 

Local 

15 

Glll.23-1.23 

111.23 

-1.24 

io.o3;°:g 

18.0+:? 

33il7 

3il5 

-2il5 

Perseus 

5 

G111.25-0.77 

111.26 

-0.77 

10.04;°;!^ 

18.1+:? 

0i7 

-lli7 

-15i9 

Perseus 

5 

NGC 7538 

111.54 

0.78 

^■^+-0.07 

14.8+:? 

19i2 

-23i3 

- 10 i 2 

Perseus 

15 

PZ Cas 

115.06 

-0.05 

V-OJ_o.i 2 

17.U_o.9 

17i3 

7i3 

-6i4 

Perseus 

13 

L 1287 

121.29 

0.66 

0 qc +0.02 

O.OJ- 0.02 

r 1 + 0.2 
' 1 - 1 - 0.2 

10i3 

-10i4 

-3i2 

Local 

21 

W3(OH) 

133.95 

1.06 

Q 7Q+0.03 
^■'^-0.03 

8 7+0.1 

19i3 

-lli 2 

li 2 

Perseus 

28 

S Per 

134.62 

- 2.20 

lO.lStHs 

Q "T + O.S 
^+-0.3 

3i3 

-14i3 

- 6 i 2 

Perseus 

3 

L 1448 

158.06 

-21.42 

0 ^o+ 0.02 

o.j:)-o.oi 

0.5 

-16i6 

-15i4 

-4i4 

Local 

10 

NGC 1333 

158.35 

-20.56 

c 54+0.02 

0.5 

-10i3 

0i2 

4i2 

Local 

8 

Hubble 4 

168.84 

-15.52 

8.46 

0.2 

-12i2 

liO.4 

-9il 

Local 

26 

IRAS 05168H-3634 

170.66 

-0.25 

10.19++ 

1 7+0.2 

1 + -0.1 

8i2 

-13i6 

-7i8 

Perseus 

23 

HP-Tau/G2 

175.73 

-16.24 

8.48 

0.08 

-18i2 

-liO.2 

Oil 

Local 

27 

T-Tau/Sb 

176.23 

-20.89 

8.47 

0.1 

-18i0.9 

12i0.1 

-2i0.4 

Local 

14 
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Table 5. (Continued) 


Source 

1 

b 

R /3t 

u 

V 

W 

Arm 

Ref.S 


deg 

deg 

kpc deg 

km 

km 

km 



G183.72-3.66 

183.72 

-3.66 

q qi +U.UJ q ^+u.ui 

^■^^-0.03 ^-^-0.01 

0±5 

1±10 

3±10 

Perseus 

5 

IRAS 06061+2151 

188.79 

1.03 

1/A 00+0.13 1 '7+0.1 

1^-7:)_0.12 

12±1 

-22±2 

-11±1 

Perseus 

18 

IRAS 06058+2138 

188.95 

0.89 

1U.U/_0.10 ^-O-O.l 

6±1 

-14±3 

4±2 

Perseus 

19 

S252 

188.95 

0.89 

IAJ^i+0.03 1 O+0.02 

10.41 _Q Q3 l-O_0.02 

-2±3 

-9±0.6 

-2±0.2 

Perseus 

20 

G192.16-3.84 

192.16 

-3.84 

Q oi +0.10 1 Q+0.1 

^■^^-0.08 ^-^-0.1 

3±3 

-8±2 

3±1 

Perseus 

25 

S255 

192.60 

-0.05 

Q RQ+0-0’^ 9 0+0.07 

^■^^-0.06 ^-^-0.08 

3±2 

3±12 

3±7 

Perseus 

21 

Orion 

209.01 

-19.39 

8.67 -1.2 

-1±4 

-1±3 

6±3 

Local 

7 

G229.57+0.15 

229.57 

0.15 

11 00+U.2c5 -i^o+O.b 

ll.a9_o.20 ^ /-^-o.T 

15±12 

-13±14 

-10±15 

Perseus 

5 

G232.62+1.00 

232.62 

1.00 

9 44+0-07 _o 1+0.4 

^■++-0.07 °-^-0.4 

2±3 

-3±3 

1±2 

Local 

20 

G236.81+1.98 

236.82 

1.98 

10.33;^-?? -14.4++ 

-7±7 

-2±7 

-3±7 

Perseus 

5 

VYCMa 

239.35 

-5.07 

8.98;h6 

0±2 

-9±3 

-3±2 

Local 

4 

G240.31+0.07 

240.32 

0.07 

11.90;°+ -22.9 tli 

-11±6 

4±7 

-14±7 

Perseus 

5 

IRAS 07427-2400 

240.32 

0.07 

1 1 Q7+0-75 oo 1 +2.1 

11-^/_0.54 ^^-^-2.5 

-8±8 

0±9 

1±11 

Perseus 

1 

DoAr21/Ophiuchus 

353.02 

16.98 

8.21 -0.1 

2±5 

-9±1 

6±2 

Local 

14 

Sl/Ophiuchus 

353.10 

16.89 

8 .22+°i -0.1 

6±5 

-2±1 

-4±2 

Local 

14 

Sun 

- 

- 

- 

C/0=11.1 + 

Vo=12.24t2 

lVo=7.25t°:® 


24 

Average-1 




-1±1 

-6±1 

1±1 

Local 


Average-2 



- „-i /r. -J C. 

8±3 

-9±2 

-4±1 

Perseus 



peculiar motions. 

^Galactocentric azimuth [deg] (see text). 

^ W49N and G48.60+0.02, located in the Perseus arm of the first Galactic quadrant, were excluded from the averaging (see text). 

^References: (l)This paper; (2)Ando et al. 2011; (3)Asaki et al. 2010; (4)Choi et al. 2008; Zhang et al. 2012a; (5)Choi et al. 2014; (6)Dzib et al. 2010; (7)Hirota 
et al. 2007; Kim et al. 2008; Menten et al. 2007; (8)Hirota et al. 2008a; (9)Hirota et al. 2008b; (lO)Hirota et al. 2011; (ll)Honma et al. 2012; (12)Imai et al. 
2012; (13)kusuno et al. 2013; (14)Loinard et al. 2008; (15)Moscadelli et al. 2009; (16)Moscadelli et al. 2011; (17)Nagayama et al. 2011; (18)Niinuma et al. 
2011; (19)Oh et al. 2010; (20)Reid et al. 2009b; (21)Rygl et al. 2010; (22)Rygl et al. 2012; (23)Sakai et al. 2012; (24)Schonrich, Binney, gc Dehnen 2010; 
(25)Shiozaki et al. 2011; (26)Ton'es et al. 2007; (27)Ton-es et al. 2009; (28)Xu et al. 2006; (29)Xu et al. 2009; (30)Xu et al. 2013; (31)Zhang et al. 2012b; 
(32)Zhang etal. 2013; 
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Fig. 5. (a)-(c) Peculiar motions {U, V, W) as a function of Galactocentric distance {R) shown in figures (a), (b) and (c), 
respectively. Red lines in (a) and (b) were fitted to the sources located in the Perseus arm of the second or third Galactic 
quadrant (see text). Red circles show the sources in the Perseus ann of the second or third Galactic quadrant, and blue ones 
display the sources in the Perseus arm of the first Galactic quadrant. Black circles represent the sources in the Local arm. 
(d) U vs. V plane. 
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4.2. Observational indication of the density-wave 
theory 

Using the obtained peculiar motions and previ¬ 
ous VLBI astrometry results, we show U and V 
of the Perseus and Local arms as a function of 
Galactocentric distance in figures 5a and 5b, re¬ 
spectively. Interestingly, linear trends might be 
seen in figures 5a and 5b for the Perseus arm, 
but except for W49N and G48.60 h- 0.02 located 
in the Perseus arm of the first Galactic quadrant. 
For W49N and G48.60 h- 0.02, we do not have a 
clear explanation for the deviations from the lin¬ 
ear trends. As noted in the previous section, more 
astrometry results for the Perseus arm in 28.8 deg 
< /3 < 81.4 deg, corresponding to 8.16 kpc < R< 
10.02 kpc (see table 5), are crucial for explaining 
the deviations. 

According to Russeil (2007), the linear trend of 
V components was previously confirmed in stel¬ 
lar rotation velocities (see fig. 7 in Russeil 2007), 
which is basically consistent with our result of the 
gaseous V components. Russeil (2007) regarded 
Galactic corotation radius as the place where V 
component is 0 km s“^, since the density-wave 
theory predicted that directions of the peculiar mo¬ 
tions {U and V) are changed inversely at inner and 
outer corotation radii (e.g., Mel’nik et al. 1999). 
Based on the assumption, Russeil (2007) deter¬ 
mined a corotation radius of 12.7 kpc by a linear 
equation fit toward the stellar V components. In 
the same way, we conduct the linear fits toward 
both U and V components, but with W49N and 
G48.60 h- 0.02 excluded from the fits. As a result, 
we obtain corotation radii of 11.1 ± 0.3 and 12.4 
± 0.5 kpc in U and V components, respectively. 
These VLBI astrometry results can be understood 
with the assumption of the density-wave theory, al¬ 
though more astrometry observations around the 
corotation radius are required to accurately deter¬ 
mine the corotation radius. 

To evaluate whether the astrometry results for 
the Perseus arm are affected by systematic errors 
(e.g., the assumed Galactic constants, peculiar so¬ 
lar motions, and Galactic rotation curve), we com¬ 
pare the peculiar motions of the Perseus arm with 
those of the Local arm as shown in table 5 and fig¬ 
ure 5. Averaged disk peculiar motions of (Lmean, 
Knean) = (8 ± 3, —9 ± 2) km s“^ are obtained 
for 25 sources in the Perseus arm (W49N and 
G48.60 h- 0.02 are excluded from the averaging). 


which are deviated from those of (t/mean, Vmean) 
= (—1 ± 1, —6 ± 1) km s“^ for 32 sources in the 
Local arm with a significance of ^ l-cr. Note that 
the obtained errors are the standard errors. The dif¬ 
ference between the Perseus and Local arms means 
that relative offsets between them in U and V com¬ 
ponents are real, although absolute values of U 
and V can be changed depending on the system¬ 
atic errors. Figure 5d also shows the relative offsets 
between the Perseus and Local arms in [/ vs. U 
plane. Most sources in the Perseus arm of the 
second Galactic quadrant are located on the lower 
right side of the figure (see also table 5), while most 
sources in the Local arm are located around the ori¬ 
gin, except for some outliers showing U ~ —20 km 
s“^. 

To explain the peculiar motions of the Perseus 
arm in detail, we will compare these with an ana¬ 
lytic gas dynamics model in the following section. 

4.3. Comparisons between observed peculiar 
motions and an analytic gas dynamics model 

Based on Pihol-Ferrer et al. (2012) and (2014), 
we compare the observed peculiar (non-circular) 
motions listed in table 5 with an analytic gas 
dynamics model in this section. Detailed explana¬ 
tions for the model and the comparison procedure 
are summarized in Appendixes 2 and 3. 

4.3.1. The model of the Milky Way 

Through the least-squares fits described in the 
Appendix 3, we determine the reduced values 
ixl) in various combinations of model parameters 
as listed in tables 6 and 7. From models ID 1 to 
ID 4, the mode m = 2 is fixed. The models ID 3 
and ID 4 (in the cases of i = 15.0° and 20.0°) show 
lower Xu values compared to ID 1 and ID 2. From 
ID 5 to ID 8, the mode m = 4 is fixed, and ID 8 (in 
the case of i = 20.0°) shows the minimum xt value 
in the four models. 

Based on all the xt values, we regard the models 
ID 4 and ID 8 as good models in this paper. We 
note that these results are first step results, since 
we have to sophisticate both the fitting procedure 
(i.e., the least-squares fit) and the analytic model. 
For instance, we assumed the constant pitch angle 
and amplitude for the spiral potential, which may 
not represent the real Galaxy (e.g., a bifurcation or 
spur of the spiral arm). However, the sophistication 
of the model is beyond the scope of this paper. 

Figure 6 shows results of the least-squares fits 
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Table 6. Summary of the asymmetric potential parameters used in this paper. 


Parameter 

Dimension 

Notes 


(kms“^) 

Amplitude of an asymmetric potential 

-^ref 

(kpc) 

Reference position of an asymmetric potential 

i 

(deg) 

Pitch angle of an asymmetric potential 

Up 

(km s“^kpc“^) 

Pattern speed of an asymmetric potential 

m 


Mode of an asymmetric potential 

A 

(km s“^kpc“^) 

Damping term for gas motion 

s 

(km s“^kpc“^) 

Corotation softening parameter 


Table 7. Results of the least-squares fits for the asymmetric potential parameters.* 


Model 

ID 

#of 

sources 


-^ref 

i 

Up 

m 

A 

£ 

X^/d.o.f 

Memo 

1 

27 

15.3 

9.46±0.01 

5.0 

24.1±0.1 

2 

27.0±2.5 

0.99±0.03 

549/50 

30% of Eo 

2 

27 

28.1 

9.39±0.02 

10.0 

17.8±0.8 

2 

24.4±1.9 

3.2±0.2 

339/50 

50% of Eq 

3 

27 

34.7 

9.10±0.04 

15.0 

17.9±0.4 

2 

18.9±1.7 

4.2±0.1 

321/50 

50% of Eg 

4 

27 

36.2 

9.16±0.05 

20.0 

10.1±0.2 

2 

6.1±0.3 

13.2±0.3 

319/50 

40% of Eq 

5 

27 

12.5 

9.74±0.01 

5.0 

24.1±0.1 

4 

7.7±0.6 

0.98±0.03 

684/50 

40% of Eq 

6 

27 

19.8 

9.22±0.01 

10.0 

10.1±0.3 

4 

1.7±0.2 

8.8±0.3 

446/50 

50% of Eg 

7 

27 

24.6 

9.01±0.02 

15.0 

11.5±0.3 

4 

5.5±0.3 

9.3±0.3 

349/50 

50% of E© 

8 

27 

25.6 

9.48±0.03 

20.0 

16.0±0.1 

4 

7.0±0.3 

7.0±0.2 

306/50 

40% of E© 


* Bold values, V i, and m, are fixed based on previous researches (see text). Other parameters are searched within realistic values 
(see also text). is searched in a range corresponding to the phase between —tt and tt (rad) of the spiral potential from i?ref = 
10.4 kpc with an increment of 0.1 kpc (see text). f!p is searched between 10 and 30 (km s“^ kpc“^) with an increment of 0.5 (km 
s“^ kpc~^). A and e are searched in the same range between 1 and 30 (km s“^ kpc~^) with an increment of 0.5 (km s“^ kpc~^). 
V^a can be converted to surface density (S). The converted values are listed in the last column (see text). 


especially for ID 4 and ID 8. Figure 7 checks 
convergences of the model parameters for ID 8. 
A (gas friction term) is not converged well in fig¬ 
ure 7. However, A is just an artificial parameter 
and not an observable parameter. Figures 8b and 
8 c show comparisons between observed peculiar 
motions (fig. 8a) and modeled U V with ID 
4 and ID 8 applied, respectively. Averaged pecu¬ 
liar motions of (C/mean, I4iean) = (8 ± 3, -9 ± 
2) km s“^ are determined for 25 sources in the 
Perseus arm in fig. 8a (W49N and G48.603-0.02 
are excluded from the averaging). Using the mod¬ 
els ID 4 and ID 8, the systematic peculiar mo¬ 
tions are canceled out as shown in figures 8b and 
8 c, respectively. In fact, averaged residuals as 
{Uohs - Uynodeh I4bs “ Knodei) are determined to 
be (-3 ± 3, 0 ± 2) km s”! and (-2 ± 2, -2 ± 2) 
km s“^ with ID 4 and ID 8 applied, respectively. 

Figures 9 and 10 display the gas orbit models 
generated from ID 4 and ID 8, respectively. 
Observed and modeled peculiar motions for the 


Perseus arm are superimposed on the gas orbit 
models (see figures 9 and 10), and also amplitudes 
of the spiral potentials (ID 4 and ID 8) are repre¬ 
sented in the figures. Interestingly, the gas orbits 
of ID 8 concentrate around the sources located in 
the Perseus arm, meaning that our model (ID 8) 
based on the observed peculiar motions can make 
a dense gas region precisely. 

4.3.2. Ojfset between gas and spiral potential 
model 

As clearly seen in fig. 10, a dense gas region 
traced by VLBI astrometry results is located down¬ 
stream of the spiral potential model, which is con¬ 
sistent with Pifiol-Ferrer et al. (2012) and (2014), 
but not consistent with the galactic shock proposed 
by Fujimoto (1968) and Roberts (1969). The differ¬ 
ence of the gas distributions may originate in a dif¬ 
ference between analytic (Pifiol-Ferrer et al. 2012 
and 2014) and numerical (Fujimoto 1968; Roberts 
1969) solutions. 
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Fig. 6. distributions for each model parameter listed in table 7. (a)-(d) The x^ distributions for model ID 4. (e)-(h) The 
distributions for model ID 8. 

Observationally, the dense gas region at the cated downstream of the spiral arm, while smaller 
downstream was confirmed in Egusa et al. (2011). clumps were located upstream of the spiral arm. 
Egusa et al. (2011) conducted CO(J=1-0) obser- In addition, the massive GMCc were traced by the 
vations with CARMA toward the nearby grand- brightest HII regions. 

design spiral galaxy M51. They found that mas- Using a dynamical model of the Milky Way (fig. 
sive Giant Molecular Clouds (GMCs) were lo- 10), we calculate the time needed to pass the off- 
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310 


320 


330 


310 


320 


Model ID 8 


10 15 20 25 

fip (km s’’ kpc’’) 

(a) 


30 



15 20 25 

fip (km s’’ kp>c‘’) 

(b) 



10 11 12 
Rref (•'PC) 


1O 15 20 

X(km s’’ kpc’’) 


Fig. 7. distributions in various combinations of the model parameters listed in table 7, especially for model ID 8. Solid, dashed 
and dotted contours represent x^ values of 310, 320 and 330, respectively. 


set between the dense gas region and the bottom 
of the spiral potential model for estimating the typ¬ 
ical time scale of the star formation. The calcu¬ 
lated time scale is 46 Myr, which is longer than 
a star-formation time-scale of < 10 Myr proposed 
by Egusa et al. (2009). Also, 46 Myr is longer than 
the free fall time (f//) for GMCs and a lifetime of 
massive GMCs, which are ^ 10 Myr (see fig. Id 


in Koyama & Inutsuka 2000 with a number den¬ 
sity of 10^ cm“^) and ^ 30 Myr (Stabler & Palla 
2005), respectively. Note that gas depletion times 
are much longer than f f / in the Galactic molecular 
clouds (e.g., see table 4 in Evans et al. 2009). 

On the other hand, 46 Myr might be explained, if 
one assumes a magnetically supported cloud where 
a timescale of magnetic flux loss Hb) is 10—100 













18 


N. Sakai et al. 


[Vol., 


60 

-1-1-1-1- 1 - 


Perseus arm in the • 

40 

2nd or 3rcl Galactic quadrant 


Perseus arm in the 1st • 


Galactic quadrant 

^ 20 

— 

1 


c n 


C U 


-20 


-40 

- 

-60 

_1_1_ 1 _1_ 1 _ 


-60 -40 -20 0 20 40 60 


U{km s'^) 



60 
40 

20 

0 

-20 
-40 
-60 

-60 -^0 -20 0 20 40 

^obs " ^ ^ 


-1-1-1 I r 

Model ID 4 




2nd or 3rd Galactic quadrant • 
1st Galactic quadrant • 

_I_I_L__I_u 



(a) 



60 

40 

20 
0 

-20 
-40 
-60 

-60 -40 -20 0 20 40 60 

^obs ~ ^nodel ^ 1 

(c) 


- 1 - 1 I- 1 r 

Model ID 8 





2rKi or 3rd Galactic quadrant • 
1st Galactic quadrant • 

_ \ _I_I_I_ 


(b) 


Fig. 8. (a) Peculiar motions plane (U vs. V) for the sources in the Perseus arm (see table 5). Red circles show the sources 
in the Perseus arm of the second or third Galactic quadrant, and blue ones display the sources in the Perseus arm of the first 
Galactic quadrant, (b) Residuals between the observed peculiar motions and the models generated from ID 4 in table 7. (c) 
Same as (b), but from ID 8. 


times longer than the free fall time {ts ^ f//) (see 
fig. 3 in Nakano et al. 2002). Observationally, Sato 
et al. (2008) determined a dynamical age of ^ 20 
Myr for a star-forming region associated with the 
NGC 281 superbubble with VLBI astrometry. 


Based on the above discussions, the longer time 
scale (46 Myr) derived from a dynamical model of 
the Milky Way might be explained, if star forma¬ 
tion is triggered downstream of the spiral potential 
model, or if molecular clouds are magnetically sup- 


60 
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X [kpc] 

gas orbit model • Local arm —<*— Other region •- 

Perseus arm ♦- Perseus arm in the 1st quadrant 


Fig. 9. The gas orbits model generated from model ID 4 listed in table 7. Filled circles show previous VLBI astrometry 
results, and different colors mean different regions (e.g.. Red = Perseus arm in the second or third Galactic quadrant; Black 
= Local arm; Blue = Perseus arm in the first Galactic quadrant; Green = others). Black and red an'ows show observed 
and modeled peculiar (non-circular) motions, respectively. Note that we assumed a flat rotation model (i.e., 0{R) = 0o), 
Galactic constants of {Ro, 0o) = (8.33 kpc, 240.0 km s“^), and peculiar solar motions of (Uq, Vq, Wq) = (11.1 it 1, 
12.24 it 2, 7.25 it 0.5) km s“^ to derive the peculiar motions. Color bar represents amplitude of the spiral potential model 
(ID 4). Green dotted circle shows the Inner Lindblad Resonance (ILR). Solar position is (X, Y) = (0, Ro) kpc. 


ported. 

4.3.3. The "revised kinematic distance” with our 
model 

Here we model loci of the Perseus arm on the 
I — V diagram of CO(J=1-0) emission (Dame et al. 
2001) using the model ID 8 in table 7. Generally, 
the observed line-of-sight velocity (li.sR) can be 
written by 


VisR = Rosinl (2) 

with the assumption of circular motion. Here 0 
is the rotation speed at a Galactocentric distance 
of R, I is Galactic longitude, and (Rq, ©q) are the 
Galactic constants. We modify eq. 2 with correc¬ 
tion terms about the peculiar motions as 


(km s’ ) 
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Fig. 10. Same as fig. 9, but using model ID 8 instead of ID 4. Inner and outer dotted circles represent the Inner Lindblad 
Resonance (ILR) and the Corotation Resonance (CR), respectively. 


Vlsr = Rosinl - Ucos4i,i3) 

where U and V are the peculiar motions, and (p = 
180.0° — (/ + /3). Note that /3 is the Galactocentric 
azimuth explained previously. For R we can use 
the logarithmic spiral models for the Perseus arm 
as determined in this paper and Sakai et al. (2012). 
For the peculiar motions (U and V) we can refer to 
the model ID 8. 

As a result, fig. 11a shows the loci for the 
Perseus arm with or without the corrections of the 


peculiar motions. Note that we assumed a flat ro¬ 
tation model (i.e., <d{R) = ©o) in both procedures. 
The astrometry results listed in table 5 are also su¬ 
perimposed on fig. 11a. Clearly, the corrections 
with ID 8 work well to trace the sources in the 
Perseus arm on the I — v diagram compared to the 
case without the corrections. The obtained loci (as 
red curves in fig. 11a) can be used to estimate a 
distance for a source located on the loci (i.e., the 
source located in the Perseus arm). 

To evaluate an accuracy of the distance measure- 


(km s ) 
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Perseus arm 
Other region 


,TD AC 07477 




Li- V 




Flat rotation model with the log spiral of Sakai+12 
Fiat rotation model with the log spiral of this paper 
Model IDS with the log spiral of Sakai+12 
Model IDS with the log spiral of this paper 


Galactic longitude / (deg) 


(a) 



(b) 


Fig. 11. (a)VLBI astrometry results as colored circles superimposed on the longitude-velocity diagram of CO(J=1-0) emis¬ 
sion (Dame et al. 2001). Different colors mean different regions (e.g., Red = Perseus arm in the second or third Galactic 
quadrant; Blue = Perseus arm in the first Galactic quadrant; Black = other regions). The horizontal axis shows Galactic 
longitude, and the vertical one represents the LSR velocity. The solid black curve is a model for tracing the Perseus arm, 
which is based on a flat rotation model (i.e., S{R) = ©o) and the log spiral model (e.g., see eq. 1) of Sakai et al. (2012). 

The dashed black curve is the same as the solid one, but based on the log spiral model of this paper. Solid and dashed red 
curves are the same as the black curves, but with peculiar motions corrected based on model ID 8 in table 7 (see text), (b) 
Accuracies of distance measurements using the models (solid black and red curves) in fig. 11(a) relative to distance mea¬ 
surements by VLBI observations (see text). Horizontal axis shows Galactic longitude, and vertical one represents a ratio of 
both distance measurements, I -^modei i black circles were derived based on the solid red and black curves in fig. 

' ^VLBI ' ® 

11(a), respectively. 

merit about for the loci (as red curves in fig. 11a), and 

we formulate the revised kinematic distance (D) i -— 

using eq. 3 as D = RqqosI zb y B? — i?Qsin 1. (4) 

Rosinl 


R={Q + V) 


14sr + U cos(j) + 0osin/ 


Note that in the final step we just used the law of 
cosines. Using eq. 4 with the corrections of the pe¬ 
culiar motions, we can determine distances toward 


LSR radial velocity (km s‘‘) 
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the sources located on the loci in fig. 11a. Ratios of 
the revised kinematic distances to VLBI astrometry 
results are shown as red circles in fig. 11b. 

On the other hand, the same procedures are 
conducted toward the “general” kinematic dis¬ 
tances in fig. 11b (as shown by black circles). 
Especially for Galactic longitude I between 94.60° 
(corresponding to AFGL 2789) and 134.62° 
(corresponding to S Per), our model (ID 8) shows 
good corrections, although around I ~ 180° it 
still has large uncertainties. The averaged ratio in 
the case of the corrected distances is 1.03 ± 0.20 
with 94.60° < I < 134.62° for 14 sources in the 
Perseus arm, while that in the case of the kinematic 
distances is 1.45 ± 0.35. Note that the errors are 
the standard deviations. 


5. Conclusion 

We showed astrometry results of IRAS 07427- 
2400. Trigonometric parallax was obtained to be 
0.185 ± 0.027 mas, corresponding to a distance of 
5.41 ^Qgg kpc, which is consistent with the previ¬ 
ous result of 0.188 ± 0.016 mas obtained by Choi 
et al. (2014) within error. Systematic proper mo¬ 
tions were determined to be (naCosS, ns) = (— 1.79 
± 0.32, 2.60 ± 0.17) mas yr“^ with six maser fea¬ 
tures ranging 62.3 < ki,sR < 70.4 km s“^, while 
Choi et al. (2014) showed (fiaCosS, ns) = (—2.43 ± 
0.02, 2.49 ± 0.09) mas yr“^ with one maser feature 
ranging 67.1 < Elsr ^ 68.0 km s“^. Averaging of 
the multiple features is essential to accurately es¬ 
timate the systematic proper motions, since a dif¬ 
ference of observed proper motions was up to ^ 2 
mas, corresponding to ~ 50 km/s at a distance of 
5.4 kpc. 

Also, we succeeded in conducting direct (quan¬ 
titative) comparisons between VLBI astrometry re¬ 
sults and an analytic gas dynamics model for test¬ 
ing the density-wave theory, while the direct com¬ 
parisons have suffered from observational limita¬ 
tions such as insufficient spatial and velocity reso¬ 
lutions since 1964, in which the density-wave the¬ 
ory was proposed by Lin & Shu (1964). We ob¬ 
tained mainly two results from the direct compar¬ 
isons, although sophistication of the model should 
be conducted (e.g., a bifurcation of the spiral arm). 

First is the different pitch angles of gas and 
(probably) stellar spiral arms. The pitch angle of 


the Perseus arm was obtained to be 11.1 ± 1.4 
deg by VLBI astrometry results in this paper, while 
that of the spiral potential model was obtained to 
be ^ 20 deg. The difference of the pitch an¬ 
gles will be distinguishable by stellar astrometry 
(e.g., Gaia astrometry). In addition, a dense gas 
region traced by VLBI astrometry results was lo¬ 
cated downstream of the spiral potential model (see 
fig. 10), which was previously confirmed in the 
nearby grand-design spiral galaxy M51 in Lgusa 
et al. (2011). The offset between the dense gas 
region and the spiral potential model will also be 
distinguished by stellar astrometry. 

In the Gaia era, a combination of gas and stel¬ 
lar astrometry will be a powerful tool to distinguish 
several dynamical models such as the density-wave 
theory and the recurrent transient spiral proposed 
by mainly numerical simulations. The model se¬ 
lection will allow us to understand the origin and 
evolution of the spiral arm, as well as those of the 
Milky Way. 
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Here, section 3.3 is explained in detail. 

For Steps 2 and 3 (coordinates conversion and proper motion to absolute motion), we refer to the 
simple matrix multiplication shown by Johnson & Soderblom (1987) to determine full-space motions of 
the source at the solar position as below: 


■ U' ' 


"^helio 

v 

= B 

k^aCosS/n 

W' 


-1 

_1 


where 

B = T A 


with 


T = 


-I-COS0O -|-sin0o 0 


—sin^NGP 0 +cosi5ngp 


-|-sin0o —COS0O 0 


0 -10 


0 0+1 


+cosi5ngp 0 +sin(5NGP 



-l-cosaNGP 

-bsinaNGP 

0 


and 


-l-sinaNGP 

-cosqjngp 

0 


0 

0 

-fl 


A = 


-l-cosa cosi5 —sina 
-bsina cosS -fcosa 
-|-sin(5 0 


—cosa sinJ 
—sina sin(5 
-|-cos5 


Here k = 4.74057 and (U', V, W) are 3D-velocity vectors at the solar position in Galactic coordinates. 
U', V, and W' are directed toward the Galactic center, the Galactic rotation, and the North Galactic pole, 
respectively. We cite Reid et al. (2009a) to set right ascension and declination of the North Galactic pole 
(q;ngp. ^NGp) in J2000 coordinates and the position angle of the North Celestial pole {9q). 

Using the above corrections, we can write the peculiar motions at the source position with matrix 
formulas below: 


U ' 


'^helio 


Ue 


0 

V 

= c 

kflaCOsS /tt 

+ D 

Vq 

— 

e{R) 

w 


k^jLs/i: 


1+0 


0 


where 

C = D B 


and 


D = 


-l-cos/3 —sin/3 0 

-l-sin/3 -l-cos/3 0 

0 0 -bl 


Here /3 is Galactocentric azimuth (/3 is defined as 0 deg toward the Sun from the Galactic Center and 
increasing clockwise). 

For the error calculation of the peculiar motions, we show the matrix formula as 


sir 


^^helio 


Cl2 

Cl3 


= E 

(^A) [<^Lcos 5 + AaCOsJ 5^/TTf] 

+ 2/iQ cosSfis A/A 

C22 

C23 

J 




C32 

C33 


Here the elements of the matrix E are the squares of the individual elements of C, i.e., for all i 

and j. Note that the final matrix can be applied at any place of the Galactic disk, while the matrix (eq. 2) 
in Johnson & Soderblom (1987) was derived with /3 = 0 deg. Clearly, both matrices are identical if /3 is 0 
deg. 
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Appendix 2. Introduction of the analytic gas dynamics model (Pinol-Ferrer et al. 2012 and 2014) 


The model proposed by Pinol-Ferrer et al. (2012) and (2014) solves the equations of motion using the 
linear epicyclic approximation (see Binney & Tremaine 2008, p.l89), which is regarded as an extension 
of previous research related to galactic dynamics (e.g., Lindblad 1927; Lindblad 1958; Sanders & Huntley 
1976; Lindblad & Lindblad 1994; Wada 1994). 

To explain the observed velocity ellipsoid in the Milky Way, Lindblad (1927) introduced the epicyclic 
description of nearly circular orbits for stars in a circularly symmetric galaxy. Lindblad (1958) added 
a perturbing potential, showing rigid-body rotation with the pattern speed (Hp), in the theory given by 
Lindblad (1927). He pointed out that there were resonances in the perturbed system. The resonances are 
now called as the Lindblad resonances. Lindblad & Lindblad (1994), as well as Wada (1994), introduced 
gas dynamical friction (A) in the epicyclic approximation to examine gas motions. Using the theory, 
they succeeded in examining the gas motions perturbed by a bar potential without the occurrence of the 
Lindblad resonances, thanks to the introduction of A. However, another resonance called the corotation 
resonance (CR) occurred at the place where Up = il{R) (as the angular speed of star and gas) in their 
papers, which forced them to research the gas motions in inner corotation radius (where ft{R) > Up) 
for avoiding the resonance: they examined the gas motions perturbed by a bar potential around galactic 
center where il{R) > ftp. 

Following the researches described above, Pinol-Ferrer et al. (2012) and (2014) introduced artificial 
damping term (e) with the gas dynamical friction A in the epicyclic approximation to examine gas motion 
in the entire disk without the all resonances (i.e., co-rotation and Lindblad resonances). The important 
point in the new model is that they formulated the model to be able to choose any asymmetric potentials 
such as bar and spiral potentials, while comparisons between observations and the previous model pro¬ 
vided by Lindblad & Lindblad (1994) and Wada (1994) have been conducted in always bar potentials 
(e.g., Lindblad et al. 1996; Sakamoto et al. 1999; Boone et al. 2007). 

Using the new model provided by Pinol-Ferrer et al. (2012) and (2014), one can obtain steady solutions 
for gas motions perturbed by an arbitrary asymmetric potential showing rigid-body rotation with flp in 
a co-rotating flame. Mathematical descriptions for the model are summarized in the Appendix of Pinol- 
Ferrer et al. (2012), and therefore here we just explain important parts of the model. First, we write 
axisymmetric and asymmetric potentials in polar coordinates (R, 9) by 

$(i?,0) = $o(i?) + ^i(i?,0), (Al) 

where subscripts 0 and 1 denote axisymmetric and asymmetric potentials, respectively. For a technical 
reason, we regard the Galactic center as the origin of the polar coordinates, and also 6 is defined to be 
90° toward the Sun and increasing counterclockwise in this paper. In the situation, 9 is related to the 
Galactocentric azimuth /?, which is 0 = 90° — f3. 

Second, the asymmetric potential can be rewritten to be 

n 

= - ^ ^„(i?)cosm(0-^?™(i?)). (A2) 

m—1 

Here m is a mode of the asymmetric potential (e.g., m = 2 for a bar potential), is an amplitude of the 
asymmetric potential, and "dmiR) shows a phase of the asymmetric potential: In a bar potential dmiR) 
is constant, and for a spiral potential here we assume the logarithmic-spiral model as explained in section 
4.1. For 'i>miR) here we assume it as constant for simplicity in the same manner as Sumi et al. (2009). 
In the case of the logarithmic-spiral model, dm {R) can be written by 

TD 

dmiR) = cotiln —— -I- 90.0°. (A3) 

Rret 

Note that i and i?ref are the pitch angle and the reference position in the logarithmic spiral model. 

Third, R and 9 in the polar coordinates are also decomposed into axisymmetric and asymmetric parts 
as 
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R = ro+^ (A4) 

and 

0 = eQ + (n-n„)t + —. (A5) 

ro 

Here ^ and rj are deviations from circular motion, t is time, and H(ro) is an angular velocity of the circular 
motion. Note that the axisymmetric part tq is different from the Galactic constant Rq. Fourth, based on 
Pihol-Ferrer et al. (2012) and (2014), the equations of motion (see also Binney & Tremaine 2008, p.l89) 
are linearized by neglecting higher order terms of ^ and p as 

^ + 2X^-2nr]- AnA^ = - (A6) 

oR 

and 

rj + 2fl^ + 2Xr] + 4AA^ ~ ~1{ 36 ' (A7) 

Provided that we are not close to the corotation (where H = Op), we can replace eq. (A5) by 

0 ~ 6»o + (O - Op)f. (A8) 

Using eq. (A8), full solutions of the linearized equations (A6 and A7) at the guiding center are written as 


n 


^ ^ [dmCOsm{9 - '&m{R)) + emSmm{9 - idmiR))] > 

m—1 

(A9) 

n 

P = ^ [ffmSinm(6» - 'dm{R)) + fmCOSm{9 - , 

(AlO) 


m—1 


where the amplitudes dm, dm, and fm are functions of Galactocentric distance as explained in the 
Appendix of Pinol-Ferrer et al. (2012). Pinol-Ferrer et al. (2012) and (2014) inserted the damping term 
e in the amplitudes to damp the corotation resonance. 

Finally, we describe simple relations between observed peculiar motions {U and V) and the derived 
solutions (^ and p) as 

i=-U (All) 

and 

fi = V. (A12) 

Using equations (All) and (A12), we can compare VLBI astrometry results with the analytic model. 
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Appendix 3. The least-squares fit 


For comparisons between the analytic model explained in Appendix 1 and the VLBI astrometry re¬ 
sults listed in table 5, we search the minimum value of the reduced chi-square (xi) in analytic model 
parameters that are varied within appropriate ranges. The chi-square (x^) can be simply written as 


X 


2 


E 


Aoj 

o-t 


(A13) 


where is a measurement with an uncertainty cTi, and Aoi is the residual between the measurement 
and an expected value. We use U and V with uncertainties from table 5 for the measurements, and the 
residuals are calculated using the measurements and the analytic model. 

Table 6 displays the analytic model parameters used for a spiral potential, while we assume a flat 
rotation model (i.e., 0(i?) = 0o) and Galactic constants of (Rq, 0o) = (8.33 kpc, 240.0 km s“^) to 
describe a Galactic axisymmetric potential. There are seven parameters of which three parameters, the 
amplitude, pitch angle, and mode of the spiral potential, are fixed based on previous research. 

For the amplitude of the spiral potential, we refer to Grosbpl et al. (2004) who showed the relation 
between pitch angles and amplitudes for spiral galaxies (see fig. 8 in Grosbpl et al. 2004). The spiral 
amplitudes ranged between ~ 0 and ^ 50 % relative to disk amplitudes in Grosbpl et al. (2004), and 
therefore we fix the amplitude between 10 and 50 % relative to the local amplitude with an increment of 
10 %. Note that the amplitude is related to surface density S (see eq. 6.30 in Binney & Tremaine 2008) 
and we assume a local surface density (S©) of 63.9 Mq pc“^ based on Mcmillan (2011). 

For the pitch angle and mode of the spiral potential, we cite Vallee (2013) who listed recent results 
of spiral arms in the Milky Way. The pitch angles between 2.3 (or 5.4) and 16.5 deg were listed with a 
mode of 2 or 4 in tables 1 and 2 of Vallee (2013). Note that the pitch angle of 2.3 deg was determined 
based on only two sources associated with the Outer arm in Reid et al. (2009a), and recently the value 
was modified to 13.8 ± 3.3 deg based on six sources associated with the Outer arm in Reid et al. (2014b). 
Thus, we fix the pitch angle between 5 and 20 deg with an increment of 5 deg. Also, we fix the mode of 
2 or 4. 

Aside from the fixed values, we may not have to fix the reference position (i?ref) of the spiral potential 
based on VLBI astrometry results (i.e., i?ret = 10.4 kpc determined in section 4.1), since the VLBI 
astrometry results may not trace the bottom of the spiral potential. Therefore, we search it with an 
increment of 0.1 kpc in a range corresponding to phase between —tt and h-tt radian of the spiral potential 
from i?ret = 10.4 kpc. Note that the search range of i?ref depends on not only the pitch angle, but also the 
mode of the spiral potential. For instance, in the case of i = 5.0° and m = 4, we search between 9.7 
and 11.1 kpc, while we search R^^i between 5.8 and 18.2 kpc in the case of i = 20.0° and m = 2. 

For the pattern speed (flp) of the spiral potential, we search it between 10 and 30 km s“^ kpc“^ with 
an increment of 0.5 km s“^ kpc~^ based on previous research (e.g.. Op = 11.5 ± 1.5 km s“^ kpc“^ in 
Gordon 1978; Op = 30 km s“^ kpc~^ in Fernandez et al. 2001). 

The other two parameters (A and e) are not observables, and we search them in a relatively wide range. 
A and e are searched in the same range between 1 and 30 km s“^ kpc”^ with an increment of 0.5 km s“^ 
kpc”'^. 

After the grid (discrete) search, we determine final values with errors using the formula provided by 
Bevington & Robinson (2002) as 


X^ = 




■c, 


(A14) 


where Uj is a single parameter in the vicinity of the minimum (at a^) of the x^ distribution, is a 
dispersion of the single parameter, and the constant 67 is a function of the uncertainties ai (as explained 
in eq. A13) and the other parameters ak for k ^ j. 
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However, we note that the error obtained from eq. 
should be discussed in the spiral potential model. 
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